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Abstract. - We predict that spin- waves in an easy-plane ferromagnet have a finite lifetime at zero 
temperature due to spontaneous decays. In zero field the damping is determined by three-magnon 
decay processes, whereas decays in the two-particle channel dominate in a transverse magnetic 
field. Explicit calculations of the magnon damping are performed in the framework of the spin- 
wave theory for the XX Z square-lattice ferromagnet with an anisotropy parameter A < 1. In 
zero magnetic field the decays occur for A* < A < 1 with A* ~ 1/7. We also discuss possibility of 
experimental observation of the predicted effect in a number of ferromagnetic insulators. 



Introduction. — Zero-point fluctuations are generally 
present in quantum antiferromagnets, whereas their role 
in ferromagnetic structures is usually considered as minor 
[TH2] . Without challenging the fundamental reason behind 
this conclusion, we present here an example of a quantum 
effect specific to anisotropic ferromagnets. Namely, we 
predict that magnons in an easy-plane ferromagnet have 
a finite lifetime at zero temperature due to spontaneous 
three-particle decays. In contrast, excitations in a simple 
two-sublattice antiferromagnet remain stable at T = [3l 
[4] and acquire nonzero decay rate only above a threshold 
magnetic field [5l-flTj. 
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The low-frequency dynamics of an isotropic ferromag- 
net is governed by conservation of the total spin S to t • The 
ferromagnetic ground state and excited states can be clas- 
sified according to S z ot , where the z-axis is chosen paral- 
lel to the spontaneous magnetisation. In particular, the 
ground state has Sf ot = NS and every magnon carries 
an intrinsic quantum number AS Z = — 1. Spin conser- 
vation translates into conservation of the total number of 
spin waves by magnon-magnon interaction. Particle non- 
conserving processes, including spontaneous decays, are 
forbidden by symmetry for the Heisenberg ferromagnet. 

In the case of an easy-plane ferromagnet the sponta- 
neous magnetisation is orthogonal to the principal z-axis, 
see fig. QJl The rotation symmetry is completely broken 
and no value of S z can be assigned either to the ground- 



Fig. 1: (Colour online) Easy-plane ferromagnet in zero (a) and 
in transverse (b) magnetic field. The self-energy diagrams pro- 
duced by two-magnon (c) and three-magnon (d) decay vertices. 



state or to low-energy excitations: spin of a spin-wave 
ceases to exist. Absence of conserved physical quanti- 
ties other than the total momentum allows various pro- 
cesses that change the magnon number. Among them, 
spontaneous decays, fig. [TJ; and [TJl, determine lifetime of 
magnetic excitations at zero temperature. The necessary 
condition for damping is conservation of energy in an ele- 
mentary decay process. In the following we use the spin- 
wave theory to study spontaneous magnon decays in an 
easy-plane ferromagnet in the two cases: with and with- 
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out transverse magnetic field. In the former case both 
two- (fig. [lb) and three-particle (fig. [TH) decay processes 
are present, whereas in the latter case only three-particle 
decays are compatible with the mirror symmetry z — > — z. 

The predicted magnon decays at T — > may be ob- 
served in easy- plane ferromagncts such as K 2 CuF4 [T2TU4) . 
CsNiF 3 [IS] and CeRh 3 B 2 [S]. The second class of physi- 
cal systems relevant to our work is the lattice boson mod- 
els fTTHI!?] . The equivalence between an XY spin-1/2 fer- 
romagnet and a system of hard-core bosons was estab- 
lished a long time ago by Matsubara and Matsuda [2TJ1I2"!"] . 
Bose-Einstein condensates of cold atoms in optical lattices 
provide an experimental realization of such bosonic sys- 
tems [22] . The analogy between field-induced transitions 
in quantum magnets and the Bose-Einstein condensation 
(BEC) of particles was explored in many experimental and 
theoretical studies, for review, see [23]- The main focus 
was so far on the critical properties of magnon BEC. The 
dynamical aspects of magnon condensation have started to 
attract attention only recently [5 HTl"ll2"4"] . An easy-plane 
ferromagnet exhibits a quantum critical point in a trans- 
verse magnetic field [HJ[23], which can be mapped onto 
the BEC of magnons (25] ■ Therefore, investigation of the 
dynamics of a ferromagnetic model also helps to clarify 
universal dynamical properties of the Bose-Einstein con- 
densates. 
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Model. — In this Letter we investigate the zero- 
temperature dynamics of the quantum XXZ ferromagnet 
on a square lattice given by the nearest-neighbour Hamil- 
tonian 



Fig. 2: (Colour online) Magnon dispersion of the spin-1/2 
XXZ square-lattice ferromagnet in zero applied field along 
symmetry directions in the Brillouin zone. Solid curves 
show the harmonic energies for three different values of the 
anisotropy constant A. Dashed curves correspond to renor- 
malized magnon dispersion. For A = 0.9 the dashed curve is 
indistinguishable from the solid line. 



is rewritten in a tilted coordinate system such that the 
new z'-axis is directed parallel to the equilibrium mag- 
netisation. Bosonisation of spin operators is performed in 
the new frame using the Holstein-Primakoff transforma- 
tion mm 
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with an easy-plane anisotropy A < 1 and an arbitrary spin 
S. The model ([TJ describes magnetic properties of layered 
ferromagnet K2CUF4, which has a small anisotropy A ~ 
0.99 and an extremely weak exchange coupling between 
square planes: J' / J < 1CT 3 [IMU- 

In zero magnetic field the ferromagnetic moment is ori- 
ented arbitrarily in the xy-plane breaking the SO(2) ro- 
tation symmetry. In a finite field the magnetisation tilts 
away from the easy-plane by an angle 9, see fig. |T|d. The 
mean-field (classical) expression for the tilting angle 9 at 
zero-temperature is 



Sf = Sf ± iSf. Square roots are subsequently expanded 
to the first order in 1 /S and terms up to the fourth order 
H~y^ j SI (1) in boson operators <Zj and at are taken into account. 

i At the harmonic level the boson Hamiltonian is diago- 

nalized by the Bogolyubov transformation 



a k = u k 6 k 



The bare magnon energy for H < H c is 



(4) 



e k = AJSyj (1 - 7k ) [1 - 7k (Acos 2 9 + sin 2 9)] (5) 

with 7 k = i(cosfc x +cosfc y ), while the Bogolyubov coef- 
ficients are given by 



A k + e i^i/2 
2e k 



/ 4 - £k y/ 2 7k 

V 2e k J | 7k | 



sin0 = — , H c = 4JS(1 — A) 

He. 



(2) 



A k = 4 JS 



1 - 2 7k(l + A cos 2 9 + sin 2 9) 



(6) 



For H > H c local magnetic moments become completely 
aligned with the applied field and the SO(2) rotation sym- 
metry is restored. The transition at H = H c provides a 
simple example of the quantum critical point j2TJ[25j[26] . 

We use the spin-wave theory to calculate the damping 
of magnetic excitations in the model ([T]) at T = 0. The 
derivation for the anisotropic ferromagnet closely follows 
a similar analysis for a quantum two-sublattice antifer- 
romagnet in external field [5J[TU]. The Hamiltonian ([1]) 
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Zero-point fluctuations in the ground state {a\a,i 
are controlled by the parameter (1 — A) cos 2 9. They are 
gradually suppressed by the transverse magnetic field and 
vanish at H > H c . 

The broken SO(2) symmetry in an easy-plane ferromag- 
net is reminiscent of the breaking of the U(l) gauge sym- 
metry in the superfluid state of a Bose gas Ac- 
cordingly, the energy spectrum of a quantum XXZ fer- 
romagnet features an acoustic branch e k ~ k instead of 
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the 'normal' ferromagnetic dispersion e k ~ k , see fig. [U 
Direct expansion of (JS|) in small k yields 



> ck 
c 



c = 2JSV1 - A cos( 
15 + cos 4^ 



cos 2 6»(l - A) 
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(7) 



where </? = arctan k y / k x . The excitation spectrum remains 
gapless up to the critical field H c . The magnon velocity 
is positive in small fields and vanishes at H = H c . While 
the critical behaviour is mostly independent of the second 
sublcading term in the expansion (J7J , the dynamical prop- 
erties crucially depend on the convexity of the acoustic 
branch [27]. Considering the anisotropy /field dependence 
of the coefficient a(ip) we can distinguish three different 
regimes: 

(i) for A > 1/4, £k is a convex function of k for k — > and 
for all in-plane directions of the momentum; 

(ii) for A* < A < 1/4 and H = 0, the curvature of ek re- 
mains positive near the diagonal ip = ir/4 in the Brillouin 
zone but becomes negative along the principal directions 
¥> = 0,7r/2; 

(iii) for A < A* = 1/7 and H = 0, has negative curva- 
ture at k — > 0, however, there is a threshold magnetic field 
H* above which a(ip) becomes positive again: 



H* 



' 1-7A 
TO^A) 



(8) 



A similar change in the convexity of the acoustic branch 
in a finite magnetic field takes place in quantum antifer- 
romagnets [5]. 

In the isotropic ferromagnet (A = 1) the expression for 
the energy of a single-magnon excitation ek = 4 JS(1 — 7k) 
is exact, whereas for A < 1 the harmonic result §5§ is 
only approximate. We have calculated the lowest-order 
Hartree-Fock correction to ([5]) in zero field determined by 
quartic magnon terms [55]. Results are shown in fig. [5] by 
dashed lines. The quantum correction is momentum de- 
pendent shifting in opposite directions the acoustic branch 
and the high-energy magnons. For A = and S = 1/2, 
the magnon velocity is increased by 11%, which is com- 
parable to the 16% enhancement for the spin-1/2 square- 
lattice Heisenberg antifcrromagnet [3]. Accordingly, the 
change in the curvature of the acoustic branch takes place 
at A* w 0.24 rather than at 1/7. For A = 0.5, the quan- 
tum renormalization does not exceed 3-4% for any k. In 
the following we disregard completely the real part of the 
quantum correction to the spectrum of the XXZ ferro- 
magnet and focus on the magnon damping given by Im e^. 

Two-magnon decays. — In the canted magnetic 
structure at < H < H c , the principal nonlinear interac- 
tion is provided by the cubic term: 



n 3 = j\ 



A)sin2^afa.(a, + at) . (9) 




Fig. 3: (Colour online) Magnon damping in the two-particle 
decay channel for the XXZ-ferromagnet with A = 0.5 in trans- 
verse magnetic field. 



After transforming to Bogolyubov bosons one obtains two 
types of cubic vertices [5] . The magnon damping at T = 
is determined by the two-magnon decay vertex expressed 
in our case as 



L=^y k (q)(6t6t_ q 6 k + h. c .) , 



2VN 



(10) 



Vk(q) 



k,q 

Hcos6 



[7k (Uk + ^k)(u q Wk-q + U q Uk-q) 



+ 7q(«q + «q)(WkWk-q + 'U k 'Uk-q) 

+ 7k-q(«k-q + t'k-q)(wk"q + Wkfq)] ■ 

The decay vertex vanishes at H = and H > H c . In 
the high-field phase this is due to the spin conservation. 
In zero field, magnons have no well-defined spin but still 
preserve the odd parity under z — > — z. The parity conser- 
vation forbids the two-particle decays in zero field, though 
the three-particle decay processes are still possible. 

In the Born approximation, the magnon decay rate in 
the two-particle channel is given by the imaginary part of 
the diagram in fig. [TJ;: 



rk = 5E y k 2 (q)^k-e q -e k - 



(ii) 



Since Vk(q) = O^S 1 / 2 ) the decay rate ([TT]) is indepen- 
dent of the spin value Tk = 0{J). 

Spontaneous two-particle decays are allowed if the en- 
ergy conservation condition 



fk 



£k- 



(12) 



is satisfied for a given magnon dispersion. For the acous- 
tic mode ([7]) solutions of eq. (fT2]) exist only for positive 
values of the coefficient a [TUll2"7] . Following the preceding 
analysis, we conclude that the low-energy magnons in the 
model ([T]) are kincmatically unstable for A > 1/7 already 
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in vanishingly small field. Expanding the vertex (jlOp in 
small momenta we explicitly obtain for the acoustic mode 



Scaling I\ oc fc 3 is characteristic for two-dimensional mod- 
els OE] . The decay rate for a 3D anisotropic ferromagnet 
behaves similar to the phonon damping in the superfluid 
4 He: T k oc fc 5 H7]. 

Away from the k — > limit, a part of the Brillouin 
zone with unstable magnons can be determined numeri- 
cally by solving (|12[) for various incoming momenta. For 
A > A* the decay region occupies a finite area already at 
H — > and spreads out quickly over the entire Brillouin 
zone with increasing magnetic field. The magnon decay 
rate for arbitrary momenta is obtained by numerical in- 
tegration of cq. (jlip . Tk along two symmetry directions 
in the Brillouin zone is shown in fig. [3] for A = 0.5 and 
several values of the transverse magnetic field. The decay 
rate is strongest at intermediate fields H ~ 0.5H C , though 
it remains significant up to H = 0.9H C . The decrease 
of Vk(q) as H — > H c is partly compensated by a grow- 
ing phase-space volume of solutions of eq. (fl2j) . Jumps 
and peaks in the momentum dependence of I\ seeing in 
fig. [3] arc respectively produced by the decay thresholds 
and the logarithmic singularities in the two-magnon den- 
sity of states of a 2D model [29] . 

A finite lifetime of low-energy spin waves in an easy- 
plane ferromagnet in applied field was briefly discussed in 
[2"5] . Overall, physics of two-particle decays in the canted 
ferromagnetic state is quite similar to the field-induced 
magnon decays in quantum antiferromagnets [5"Ml0j. In 
particular, the asymptotic expression for I\ (|13|) remains 
valid both in the present case and for the square-lattice 
Heisenberg antiferromagnet. Also, the self-consistent 
treatment of the decay processes beyond the Born approx- 
imation removes the 2D Van Hove singularities in I\ [10] . 
Away from peaks, the decay rate is not significantly mod- 
ified by higher-order processes and the Born approxima- 
tion results presented in fig. |3] may serve as a guidance for 
future experimental tests. 

Three-magnon decays. — We now turn to zero-field 
magnon decays, which distinguish anisotropic ferromag- 
nets from collinear antiferromagnets. At H = the two- 
magnon decays (fT0|) arc forbidden by symmetry. Still, for 
A > A* the energy conservation condition (|12p can be sat- 
isfied in a certain region in the momentum space. An ex- 
ample of the kinematic threshold for two-particle decays 
is shown in fig. |4] for A = 0.5 and H = 0. As was dis- 
cussed by Harris and co-workers [3] , if magnons are stable 
against two-particle decays in the entire Brillouin zone, 
then three-particle decay processes 

£k = e P + e q + £k- P -q (14) 

are also energetically forbidden. Here, we encounter an 
opposite situation: conservation of energy and momentum 
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Fig. 4: (Colour online) Three-magnon decay region (shaded 
area) and kinematic threshold for two-particle decays (dot- 
dashed line) in the first quadrant of the Brillouin zone for the 
XXZ ferromagnet (A = 0.5) in zero applied field. 

for the two-magnon processes (fl"2"|) implies that decays in 
the three-particle channel (fT4"]) are possible as well [50] . 
To our knowledge, quantum three-particle decays have not 
been systematically studied in the literature. Therefore, 
apart from specific model results we will also discuss below 
a few general aspects of three-particle decay processes. 

The three-magnon decay region in the case of A = 0.5 
and if = is shown in fig.|4]by shaded area. Its boundary 
can be determined in a way similar to the two-particle de- 
cay threshold boundary [27l[29]. In particular, velocities 
of three created quasiparticles are equal to each other at 
the decay threshold. In fact, a stronger condition of equal 
momenta for the decay products holds for the model dis- 
persion (J5]): 

£k = 3e k / 3 . (15) 

The three-particle decays typically occur in a larger part of 
the Brillouin zone, which encloses the two-magnon decay 
boundary. 

The interaction part responsible for three-particle de- 
cays is given by 

V* = £ ^(p,q)(WL P -A + h - c -) ' ( 16 ) 

k,p,q 

Derivation of the vertex function £/k(p 7 q) follows the same 
route as the preceding analysis of the cubic terms. We 
skip, therefore, further technical details and avoid present- 
ing a cumbersome expression for [7k(p, q). The magnon 
damping is given by the imaginary part of the diagram in 
fig.HH: 

r k = g ^(P, q) 5 ( £ k - e p - e q - e k _ p _ q ) . (17) 
p,q 

Since £7k(p, q) = O(J), the decay rate (fT7)l scales as I\ ~ 
J/S representing a higher-order quantum effect compared 
to the two-magnon damping (jllj) . 
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Fig. 5: (Colour online) Zero-field magnon decay rate for the 
XXZ spin-1/2 ferromagnet along the F-M line in the Brillouin 
zone. The inset: the low-energy asymptote of Fk on a double- 
log scale; the analytical result (full line) versus numerical data 
for A = 0.5 (squares). 



The momentum dependence of Tk along the diagonal in 
the Brillouin zone is shown on the main panel of fig. [5] for 
three different values of the anisotropy parameter A. Be- 
cause of the spin dependence of eq. (fTT)). we used 5=1/2 
for illustration. For other values of spin Tk is further re- 
duced by a factor 1/(25). The decay rate behaves linearly 
Tk ~ Afc in the vicinity of the decay threshold boundary, 
whereas the Van Hove singularity inside the continuum 
yields a nonanalytic contribution Tk ~ Afc In |Afc|, which 
shows up as a shoulder-type feature on all curves Tk in 
fig. [5] For both types of anomalies one finds an extra 
factor Afc compared to the behaviour of the two-particle 
decay rate [29] . 

At small momenta the interaction vertex (|16j) exhibits 
a complicated non-analytic behaviour, similar to the mo- 
mentum dependence of spin- wave interactions in quantum 
antiferromagnets (3j[3TJ . In particular, Z7k(p, q) can have 
a finite limiting value depending on the order in which 
k, p and q are taken to 0. Nevertheless, for the on-shcll 
processes (|14j) the vertex acquires a normal hydrodynamic 
form [27] 



J7k(p,q) oc \Jk-pqk' , k' = |k - p - q| 



(18) 



and vanishes as k — > 0. Substituting (|T5|) into cq. ([T7|) 
we obtain the power-law asymptote for the decay rate at 
low energies: Tk oc fc 7 . The derived power-law behaviour 
is checked on the inset of fig. [5] The deviation between 
the straight line and the data points at small k is due to 
the representation of the delta-function by a finite-width 
Lorcntzian. In the three-dimensional case the power-law 
exponent changes to n = 11. Both cases are included in 
the general expression 
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Fig. 6: (Colour online) Zero-field magnon decay rate at the M- 
point [k = (n, ir)] in the Brillouin zone versus the anisotropy 
parameter A. 



This is again to be compared with a magnon damping in 
the two-particle channel: Tk ~ k 20 ^ 1 [61129] . 

An interesting aspect of the obtained results is the de- 
pendence of Tk on anisotropy. At small momenta the 
damping decreases with an increase of A following the sup- 
pression of the decay vertex f/k(p, q). In the vicinity of 
the Brillouin zone boundary the tendency is reversed and 
Tk is primarily determined by the expanding decay region 
and a growing phase-space volume of the decay solutions 
of eq. ([F4]) . Eventually, decays disappear completely at 
A = 1. The details of such a non- monotonous dependence 
are shown in fig.|6]for the magnon damping at the M-point. 
Summarizing the above results we conclude that the easy- 
plane anisotropy in the range A = 0.5-0.85 is most favor- 
able for an experimental observation of the three-particle 
decays. 

In a finite magnetic field (fig. QJ>) the decay rate is 
given by a sum of contributions from different decay chan- 
nels. Therefore, for not too small H, the magnon life- 
time is dominated by two-particle decays. Still, three- 
magnon decay processes, being less restrictive on the value 
of the initial momentum, will promote decays beyond the 
two-particle decay region. Higher-order two-magnon pro- 
cesses, such as a self-energy insertion in the inner lines 
of the diagram in fig. [Ifc, can be regarded as an effective 
renormalization of the three-particle decay vertex V4 and 
have a similar effect on the promotion of magnon decays 
|10j . As a result, magnons in the entire Brillouin zone ac- 
quire finite lifetimes. Note, that such a behaviour is not 
entirely universal. For the triangular antiferromagnet two- 
and three-magnon decay regions coincide and the magnon 
damping occurs only in a part of the Brillouin zone [29] , 

Conclusions. — We have theoretically investigated 
zero-temperature magnon decays in an easy-plane ferro- 
magnet with the exchange anisotropy. The phenomenon 
of spontaneous magnon decays depends only on symmetry 
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(absence of spin conservation) and kinematics (curvature 
of the spectrum). Hence, decays should be also present 
in ferromagnets with the single-ion anisotropy and differ- 
ent lattice structures. The most significant damping is 
produced by two-particle decays, which are induced by a 
transverse magnetic field H < H c . In regard to a possible 
experimental observation of the predicted effects, square- 
lattice ferromagnet K 2 CuF4 with A « 0.99 is perhaps 
too isotropic to have a sizable magnon damping. Quasi- 
one-dimensional ferromagnet CSN1F3 |15j has a moderate 
single-ion anisotropy and provides a good experimental 
system to test the predicted effects. 

Spontaneous three-particle decays is the only damp- 
ing mechanism in collinear magnetic structures. For the 
model case studied here the corresponding decay rate is 
somewhat small r k /J < 0.07 (S = 1/2). However, 
the significance of three-particle decays extends beyond 
the present model. In a collinear antiferromagnet, three- 
magnon decays arc usually impossible due to the kine- 
matic constraint [3]. However, sizable frustration may sig- 
nificantly modify the excitation energy e k . In particular, 
if an acoustic branch has an upward curvature at H = 0, 
then the energy conservation will allow three-particle de- 
cays. This scenario is realized in the vicinity of the quan- 
tum Lifshitz point. A typical zero-temperature transition 
in frustrated spin models is a second-order transition be- 
tween commensurate and incommensurate antiferromag- 
netic structures. At the transition point, the magnon ve- 
locity vanishes at least in one direction in the momentum 
space producing a positive curvature of e^. Hence, three- 
particle decays are present in the commensurate state, 
whereas two-particle decays appear in the incommensu- 
rate phase [25] . 
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